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SUMMARY 


The recently developed Johnson and Tevaarwerk fluid rheology model was used to 
investigate the traction behavior for typical traction drive contacts. Fluid shear mod- 
ulus and limiting fluid shear strength were allowed to vary over the contact area in ac- 
cordance with observed traction behavior. The influence of aspect ratio of the contact 
and the invariably present spin was investigated. Graphical solutions of the analysis 
are provided which allow for the design optimization of traction drives. Comparisons 
were made with the commonly used rigid-plastic analysis of Wemitz. 

The influence of low spin values on traction was found to be negligible; therefore, 
the simple slip analysis with the elastic- plastic model is sufficient. At moderate val- 
ues of spin the complete spin analysis with the elastic- plastic model is needed to give 
accurate traction predictions. At sufficiently high spin the Johnson and Tevaarwerk 
model may be simplified to the rigid-plastic model, and the elastic effects in the fluid 
may be neglected. Contacts with a low aspect ratio predict a superior performance in 
that they show less slip for the same degree of traction. Sideways forces due to the 
spin on the contact can be substantial and may reach 75 percent of the transmitted trac- 
tion force. Also the effect due to spin is always deleterious to the traction efficiency. 


INTRODUCTION 

Typically the power transmission in continuously variable speed traction drives 
occurs in a small, elastically deformed region on the two elements. Due to the rolling 
motion and in the presence of a fluid, a thin layer of lubricant is drawn into this contact 
area, thus separating the two bodies. 

Relative velocities between the two bodies in contact shears the lubricant. Its re- 
sistance to shear permits a transmission of force between the two bodies. The amount 
of force that can be transmitted for a given amount of slip has been the subject of many 
past and present investigations, many of them resulting in fluid models (refs. 1 to 5). 

The fluid in the gap is subjected to very high pressure gradients and pressures for 
very short times. The common behavior to high pressures is an increase in the vis- 
cous shear resistance of the lubricant (refs. 6 and 7), however, due to the extremely 
short duration of the pressure, the situation is not at all clearly understood. This is 
because there are no simple experiments that will directly measure the shear properties 



of the lubricants. It is known, however, that there are definite time delays in the re- 
sponse of these lubricant properties to sudden pressure changes (ref. 8), and it is 
thought that the lubricant behaves in a more elastic fashion at low deformation rates 
(ref. 5). 

Due to the uncertainties in the shear response of the fluid, the past design tech- 
niques for analyzing traction drive contacts were based upon the rigid- plastic concept 
as used by Wemitz (ref. 9) and Magi (ref. 10). More recently a refined traction model 
was proposed and experimentally validated by Johnson and Tevaarwerk (ref. 5) which 
considers the lubricant behavior to be viscoelastic reducing to an elastic- plastic model 
at high pressures and shear strains. 

The objective of this investigation is to (1) examine the traction behavior of a lub- 
ricant for a typical traction drive contact using the Johnson and Tevaarwerk model; 

(2) generate design charts for traction drive contacts using this model; (3) make a 
comparison between this new model and the Wemitz- Magi approach; and (4) establish 
design criteria for the application of either of these models. 
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SYMBOLS 

contact area, m^ 

Hertzian contact ellipse parameters; a is the semiaxis in the rolling di- 
rection; b is the semiaxis perpendicular to the rolling direction, m 

speed pole location, m 

rate of strain tensor, sec - * 

traction force in x- and y- direction, respectively, N 
dissipative strain rate function, sec - * 
dimensionless strain rate function 

o 

elastic shear modulus, N/m 

o 

average elastic shear modulus over contact area, N/m 

9 

compliance corrected shear modulus, N/m -1 

9 

elastic shear modulus of disk material, N/m“ 
film thickness, m 

dimensionless slip in rolling direction 
dimensionless spin 

dimensionless traction in rolling direction 
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J 5 dimensionless traction perpendicular to rolling direction 

J 6 dimensionless spin torque 

dimensionless power loss 

k aspect ratio of Hertzian contact, b/a, dimensionless 

m initial slope of experimental traction curves, dimensionless 

m' initial slope of traction curve for dry rolling bodies, dimensionless 

N normal force on contact, N 

9 

P Hertzian pressure in contact, N/rrT 

9 

P Q maximum Hertzian pressure in contact, N/irT 

T spin torque on contact, N- m/rad 

q aspect ratio dependent constant, dimensionless 

U* velocities in x-direction on upper and lower body, m/sec 

AU rolling velocity difference (slip), m/sec 

V* velocities in y-direction on upper and lower body, m/sec 

AV velocity difference in y-direction, m/sec 

AX dimensionless spin pole offset 

X, Y dimensionless Cartesian coordinates 

x, y, z Cartesian coordinates, m 

Z dimensionless stress 

<5 X spin pole offset, m 

e dimensionless speed pole location 

p maximum traction coefficient from experimental traction curve, dimension- 

less 

9 

r shear stress, N/m 

9 

t reference stress for hyperbolic sinh, N/m“ 

V 

2 

r 0 limiting shear strength of fluid, N/m 

t q average limiting shear strength of fluid over contact area, N/m 

<p equivalent strain rate parameter 

go angular velocity difference between contacting bodies normal to contact 

area, rad/sec 
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Subscripts: 


e equivalent 

ij tensor coordinates 

x, y, z Cartesian coordinates 


THEORETICAL ANALYSIS 

Typically the power transmission in variable speed traction drives occurs in highly 
stressed contact areas that operate in the elastohydrodynamic lubrication (EHL) re- 
gime. The traction characteristics of these drives therefore depend to a large extent 
on the fluid and material properties at the EHL conditions that prevail. At these con- 
ditions of pressure (=d.O^ Pa), transit time, (10 - ^ sec), and temperature (30° -> 150° C), 
Johnson and Tevaarwerk, in reference 5, found that the shear behavior of the lubri- 
cants examined was best described by the following constitutive equation: 



dt 



( 1 ) 


where 

T-j stress tensor 

e.j rate of strain tensor 

T e equivalent stress, 

V^ij~ V/2 

F(t_) nonlinear viscous function to be specified later 
G elastic shear modulus of fluid 

Equation (1) was formulated by the simple addition of elastic and viscous compli- 
ances to account for the total compliance of the fluid film. Its correctness in formula- 
tion was extensively tested and verified in reference 11 using a point- contact- type disk 
machine that permits the study of traction under both slip, spin, and slip and spin. 

To obtain the tractive capacity for a typical traction drive contact, equation (1) has 
to be solved as a series of coupled differential equations with the proper use of the 
geometric boundary conditions and the applied rate of strain tensor. 

A typical Hertzian contact as it occurs in many forms of traction drives is shown 
in figure 1 together with the nomenclature that is commonly used. For this simple 
geometry equation (1) reduces to the following: 
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d(r x /G) 

dt 


T„ • ^ 

+ ^F(r e ) = 2e x 

T e 


> 


d(r y /G) 

dt 


+ -^F(r e ) = 2e 
T e 


J 


( 2 ) 


and 



The system of shear stresses acting on the surface gives rise to tractive forces F y 
and F x . 

Figure 2 shows the side view of the contact; the following assumptions have been 
made: 

(1) Film thickness h is constant over the contact area. 

(2) Pressure distribution is Hertzian and the displacement of the effective pressure 
center from the contact center is indicated by <5 X . 

Tn most f OTITIS of tT 3 Ction driYGS, not 3,11 of the Velocities GCCUT 2 ,t OGCG Gild in TTlGSt 
cases AV = V = 0 and only slip in the rolling direction and spin are present. 

The simplified contact is shown in figure 3 together with the shear strain rate dis- 
tribution in the y-direction (fig. 3(b)) and in the x-direction (fig. 3(c)). From this fig- 
ure it is clear that 
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The system of equations (2) may now be written as 


d(T x /G) + Ix f (t ) = AIJ- coy 
dt r e 6 h 


d(r y /G) 

dt 



w (x + 5 X ) 
h 




> 


j 


( 3 ) 
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We may use the following substitutions to nondimensionalize the set of equations (3): 


X = x/a 


Z x t x /7 o 


Y = y/b 
z y = t / f o 


AX = 6/a 
Z e = r e /7 o 


where t 0 is some reference stress to be defined later. From the kinematics it fol- 
lows that x = Ut. Now, because U is constant, we obtain dt = (a/U) dX. Equation (3) 
may now be written in a simple nondimensional form as 


dz x (f n /G) Z x 

— — +— F(Z • 7 ) - = 

dX Z a e o 

e 


_a = _a /usf ab \ / AU Yy/k\ 
U U \ h /VcoVab / 


dZ v<V G > . Z v 

dX 


FiZ — 1 — — /ar\/ab\ (X AX) 

e 0 U U V h / Vk 


(4) 


V o p 

Z“. + Z~ and k = b/a is the aspect ratio of the contact. The boundary 
x y 


conditions for this set of equations are 


Z x=V° 


at X 


-f: 


under the assumption that the shear stresses are zero in the inlet region. 


(5) 


NONLINEAR FUNCTION F(t ) 

With the knowledge of the nonlinear function F(t ), the set of equations (4) may be 
solved for the stresses and hence for the traction forces. 

Various forms of the function F(r ) are possible. In reference 11 it was shown 
that a fairly general form of it is given by the hyperbolic sinh, namely, 


F(r ) = — sinh 
V 


( 6 ) 


where 

r) linear Newtonian viscosity 
r reference shear stress 

V 
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For fluids under high pressure, such as those that occur in traction drive contacts, 
it is acceptable to use the limiting case of equation (6), also known as the limiting 
shear stress model. In mathematical form the function may be written as 

F(r e ) = 0 for r e < t 0 


and 


• e„) 

F(r e )= il ^ for r g = r Q (7) 

T o 

where t 0 is the local limiting shear strength of the fluid. 

The model was suggested in reference 1 and has been used extensively in the 
analysis of traction drive contacts. The advantage in using this form of the nonlinear 
function is that we now only have one disposable parameter. When this form of the 
limiting shear stress formulation is used, equations (1) reduce to the well known 
Prandtl-Reuss equations for elastic, perfectly plastic material. 

Curves B and C in figure 4 (reproduced from ref. 11), compare actual and predicted 
traction behavior of a circular contact under combined spin and slip and using the 
elastic- plastic form of equation (1). The two parameters G and t q were obtained 
from the measured initial slope and peak traction coefficient of the zero spin traction 
curve A in figure 4. These were, combined with the kinematic information on the 
strain distribution, used to predict the combined spin and slip curves B and C. The 
results show that the elastic- plastic form of equation (1) indeed predicts the observed 
traction correctly. 

VARIATION OF THE FLUID PROPERTIES WITH PRESSURE 

In order to solve equations (4) we have to know something about the effect of pres- 
sure on the magnitude of both the shear modulus G and the limiting shear stress t q . 

A simple approach would be to ignore any variation of the G and t q due to pressure 
and to set G = G and t q = t q where the bar denotes average properties; however, it 
was shown in reference 5 that to a fair approximation both G and t q are linear func- 
tions of pressure and, therefore a better description of the two parameters over the 
contact area would be to let 
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and 


3 S, ^2 3 

t o=- t oV 1 - x - Y 


G = — G W1 - X 2 - Y 2 (8) 

2 T 

This latter type of variation of G and r Q does not necessarily give the best predictions 
of traction as is shown in figure 5. Here both types of variations of G and t q are 
shown. It may be observed that for small values of spin, the traction is predicted bet- 
ter by the averaged properties t q and G over the contact area, while at the larger 
values of spin the property variation as shown in equation (8) gives better results. 

Also, it should be noted that the small strain results are predominantly elastic and 
hence governed by G, while the large strain results are mostly governed by the plastic 
properties t q of the liquid. Therefore the following variations of G and r Q should be 
used if a better fit is wanted: 


G = G 

t = — 7 - X 2 - Y 2 

o 2 0 f 

The main reason for the apparent contradiction in the experimental variation of G with 
pressure and the variation giving the best fit is thought to be caused by disk compliance 
which is discussed in appendix A. When the present analysis is used to predict the per- 
formance of traction drive contacts consisting of the lubricant and the disks, then the 
appropriate value of G to use is that for the system, not just the fluid. This system 
modulus G can be taken as the shear modulus of the lubricant to give nearly correct 
predictions. 

The shear modulus G and the limiting shear stress t q are parameters that are 
not readily available but can easily be obtained from experimental traction curves. The 
measured initial traction slope m and the peak traction coefficient p (see fig. 6) from 
a pure slip traction curve can be used to calculate G and the limiting shear stress r Q ; 
thus, 
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( 10 ) 


G = - m 
8 



= ii2L 

xrab 


It is important that these two parameters are obtained under identical operating condi 
tions to those prevailing in the traction drive contacts, namely contact pressure, 
temperature and area, film thickness and speed. It is possible, however, to use data 
that was obtained at the same pressure, temperature and speed but at different contact 
area, aspect ratio and film thickness through the use of the analysis as outlined in 
appendix B. 

For ease of notation, various nondimensional groups are introduced now. The non- 
dimensional groups are expressed in both the fluid properties G and r D and the trac- 
tion curve parameters m and p. They are 


G>/ab AU 
U 


1 T o h 
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(dimensionless slip) 


( 11 ) 


and 


T _ G-\/ab cv x/ab 

do • 

7 o h U 


— 7 t ( — \ . ^/k (dimensionless spin) 

8 W V U / 


Equations (4) and (7) may now be written as 


( 12 ) 
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Equation (13) may be solved for Z x and Z y using various values of J^, J 3 , k, and 
AX. In the reported calculations AX has been kept zero. 


RIGID- PLASTIC MODEL 

A particularly interesting case exists for equation (13) when the shear modulus 
goes to infinity and the constitutive equation reduces to that for a rigid-plastic model. 
The equations now represent the von Mises criteria for plastic flow. This model has 
been used quite frequently in traction drive analyses for example by Wernitz in refer- 
ence 9 and more recently by Magi in reference 10. The resulting stress equations for 
this model are 
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( 14 ) 


Z - 3 V 1 - X 2 - Y 2 (e - Y Vk ) ^ 

X 2 <p 

> 

Z -3 Z± - X 2 - Y 2 (X + AX) 

^ 2 (p \/k J 

where 

J (t . Y Vi0 2 + g W 

▼ k 


and 



The parameter e has a special meaning in traction drive analysis and is often referred 
to as the speed pole location. Its meaning is made clear by looking at figure 3(c) where 
e is given by e = AU/co. If we now divide this by Vab, we obtain e/vab = e = 
Au/wVab- So the speed pole location parameter e is nothing more than the ratio of 
slip to spin in the contact. This speed pole parameter is identical to that used by Magi 
in reference 10. From the definition of and Jg we find that e = Jj/Jg. Strictly 
speaking the parameter and Jg should not be used for the rigid-plastic model be- 
cause both contain the shear modulus, however, the ratio J^/Jg forms a convenient 
independent parameter that is consistent with the work of previous researchers. Equa- 
tion (14) then, is a limiting form of equation (13) when the elastic effects are negligible 
with J-j/Jg as one of the independent parameters. 

TRACTION CALCULATIONS 

Of real interest in traction drive design are the traction, the slip, and the total 
losses in the system. Also of concern is the degree of sophistication that should be used 
to obtain correct answers. The level of the analysis is dependent upon where the elastic 
properties of the fluid are important and when the plastic properties are more dominant. 
These questions may be answered by solving equation (13) for the traction parameters. 
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In terms of the nondimen sional parameters we have for the traction forces 


l T x dA=abT o I 


Z dX dY 


= 7rabr 0 J 4 


where 


>( 15 ) 


1 

7T 



Z x dX dY = nondi men sional traction in the rolling direction 



7rabr 0 p N 


J 


Similarly, 


Jr- = ^ = nondi mensional. traction perpendicular to the rolling direction. 

° 7rahf o pN 


For the torque on the spinning contact we have 


T = 



(Ty ' x “ T x ' y)dA 


or 


T _ T j 

T 0 ab7r • Vab pN * Vab 6 


(16) 


where 


J 6 




X - Z • Y • k)dX dY 

X 
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The power dissipated in the contact can now be calculated in terms of these nondimen- 
sional parameters; namely, 

J? = J 6 * J 3 + J 4 • 1 


or in terms of real power 


J 7 = 


G 


TrVabf^hU 


(T . oo + F • AU) 


y 


(17) 


or 


j 7 =^- Vk- 

‘ 8U 


jn_ ( T 

2, T 

p N 


00 + F 

X 


AU) 


J 


where the terms in parentheses are the individual components of the power loss due to 
spin and longitudinal slip. 

The important information from the calculations to present are traction J^, Jg, 
and torque Jg as functions of slip and spin Jg. Any other quantity may be calcu- 
lated from these five parameters. As an independent variable we may select slip J^, 
based upon the elastic- plastic argument, or the speed pole location Jj/Jg as found in 
the rigid-plastic analysis. For the traction in the rolling direction it was decided 
to present the data using both slip and speed pole locations J^/ Jg as the inde- 
pendent variable since both types of graphs reveal significant information. 

The zero spin traction curves for an elastic- plastic material at various aspect 
ratios are shown in figure 7. We observe that at any slip we obtain more traction 
with the lower aspect ratio contacts. The explanation of this behavior lies in the 
fact that the mean shear strain is directly proportional to the contact dimension in the 
X- direction. The lower aspect ratio contacts have their semimajor axis in the running 
direction and therefore see a larger mean shear strain and hence the mean shear stress 
is larger. The effect shows up in the traction. 

It may be expected that, at high local strain rates, the elastic effects are not 
significant and that the rigid-plastic model may be used to solve for the tractions. Fig- 
ure 8 shows the traction for the rigid-plastic model at various aspect ratios as a 
function of the speed pole location. Again lower aspect ratio contacts appear to have a 
slight edge on traction for the same speed pole parameter. 
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To observe the effect of spin Jg on the traction J 4 for the elastic- plastic model, 
the traction J 4 was calculated with various values of spin Jg present and the results 
are shown in figure 9, These results are also presented with the speed pole parameter 
as the independent parameter in figure 10. A simple comparison with figures 7 and 8 
therefore will show the influence of both elastic effects and spin on the traction J 4 . In 
each case the connected lines represent constant values of spin Jg. 

Figure 11 shows the resulting side force traction J g for k = 1 when a contact is 
under spin. This force is the result of elastic behavior of the fluid under spin and slip. 
At low slip < 1. 0, the side force J g starts at zero for zero spin, reaches a maxi- 
mum and then decreases for increasing spin. When more slip J 4 is introduced, the 
behavior is altered somewhat and the maximum value diminishes. Figure 12 shows the 
variation of the maximum value of the side force J g as a function of the aspect ratio k. 
Generally when k is small, there are larger side forces. As can be seen from fig- 
ure 12 the side force J g can be an appreciable fraction of the traction J 4 . 

Figure 13 shows the values of nondimensional torque J g as a function of speed pole 
locations of J^/Jg. Solid lines again indicate constant values of spin Jg. The resulting 
traction curves of J 4 versus slip J-j^ (fig. 9) or speed pole location Jj/Jg (fig. 10) 
show that in both cases the traction J 4 tends toward a limiting traction curve and 
appears to be independent of spin Jg. In figure 9 the lower spin value traction curves 
are indistinguishable from the zero spin traction curve. This suggests that below a 
certain value of spin Jg the resulting value of traction J 4 is insensitive to spin and 
that the simple elastic- plastic model without spin predicts the correct tractions. From 
figure 10 it may be observed that for high values of spin Jg the curves are indistin- 
guishable from the traction curve for a rigid- plastic- like material. The suggestion here 
is that for values of spin Jg above a certain threshold value the traction results J 4 
are in fact independent of any elastic effects in the sheared material. This form of 
behavior with spin Jg could have been predicted from the nature of the constitutive 
equation (13). For small values of spin Jg the behavior would be almost elastic- like 
while for large values of spin Jg the differential equations (13) become singular in 
nature. From all of the results in figures 9 and 10 it may be observed that the two 
threshold values are dependent on the ellipticity of the contact. In traction drives we 
are generally only interested in the traction range that lies somewhere near 75 percent 
of the maximum obtainable traction. For this traction value the various threshold val- 
ues of spin Jg as a function of the aspect ratio are shown in figure 14. The three re- 
gions of influence are indicated on figure 14. They are 

Region I - elastic-plastic analysis without spin is sufficient 

Region II - elastic- plastic analysis with spin is required 

Region III - rigid-plastic analysis with spin is sufficient 
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The lines dividing the regions are the values of spin Jg corresponding to a 5 percent 
increase in the slip traction level, compared to the master traction curves in figures 7 
and 8. 


LOSS FACTORS 

Losses in the contacts of a traction drive are made up of two components (slip and 
spin losses) as was indicated in equation (17): 

Jy — (Jg X Jg + X J^) 


From the calculated results the power loss Jy can be calculated from the individual 
components and we can establish curves of Jy as a function of slip or speed pole 
location J^/Jg. An example of Jy versus J^ is shown in figure 15 for k = 1. 

It is, however, more convenient to define a loss factor LF which has a more di- 
rect meaning to us as follows: 


(LF) = 



U7r Vab r^h 


(T • u> + F • AU) 


7rabr 


G Vab / power loss \ 
T Q h \ power input/ 


= !iVk- 

8 n 


( power loss 
power input 


(18) 


The loss factor may be thought of as an equivalent creep in the power transmission di- 
rection. It is interesting to note that, when the drive is idling (i. e. , no net power out- 
put), then the loss factor is simply given by 


(LF ) 0 



(19) 
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or in terms of traction curve parameters 


(LF) =-Vk- (20) 

8 n 

and we see that the idling losses are directly related to the traction curve parameters. 
Figure 16 shows the calculated loss factor curves as a function of the dimensionless 
transmitted traction J 4 . This way of presenting the data is most useful since most 
forms of the traction drives employ a mechanism whereby the ratio of the normal load 
to the tractive force transmitted is maintained constant at values of 60 to 80 percent of 
the maximum obtainable traction coefficient. This feature means that these types of 
drives operate at essentially a constant value of J 4 much of the time, that is, along a 
vertical' line at a constant value of J 4 in figure 16. 

From the results in figure 16 it may be observed that the region of 60 to 80 percent 
of the maximum available traction, or J 4 = 0. 6 to 0. 8, is the most efficient when sub- 
stantial spin is present, though this depends somewhat on the aspect ratio k. Also in 
this same region it may be observed that the loss factor is not affected by small values 
of spin and that it is the lowest for the contacts of low aspect ratio. The insensitivity 
to small values of spin indicates that the losses in the contact are predominantly due to 
creep for these small values of spin. 

In figure 17 the loss factor at traction J 4 = 0. 75 is plotted against the spin Jg for 
the various aspect ratios considered in this work. From this figure it may be observed 
that the imposed spin has little or no influence on the loss factor at low values of spin. 

From a practical point of view it can be concluded from figure 17 that a traction 
drive may be operated with values of spin Jg up to approximately one without increasing 
the power losses in the contact. 

Figure 18 shows the loss factor plotted as a function of J^/Jg at constant values 
of spin Jg for k = 1. This figure indicates that the minimum value of the loss factor 
occurs at approximately constant values of the speed pole location e or J^/Jg in- 
dependent of the amount of spin present. This may be useful information in the design 
of angular contact bearings to minimize their losses. 


CONCLUSIONS 

The elastic-plastic approximation of the Johnson and Tevaarwerk fluid traction 
model was used to analyze the traction behavior of typical traction drive contacts. The 
two required fluid properties, shear modulus and limiting fluid shear strength, were 
allowed to vary over the contact area in accordance with observed traction behavior. 
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This investigation was particularly aimed at the role that both contact configuration and 
spin play in longitudinal traction. Graphical solutions of the analysis are provided which 
allow for the design optimization of traction drive contacts. Comparisons were made 
between the full Johnson and Tevaarwerk solutions and those which can be arrived at by 
using the commonly used rigid-plastic analysis due to Wemitz. 

The most important conclusions drawn from the analysis are the following: 

1. Under conditions of slip and spin several simplifications in the analysis may be 
used; they are (a) The influence of low spin values is negligible on the traction. Simple 
slip analysis with the elastic-plastic model is sufficient, (b) At moderate values of spin 
the complete analysis of spin with the elastic- plastic model is needed to give accurate 
traction predictions, (c) At high values of spin the elastic effects in the fluid may be 
neglected and the simple rigid-plastic analysis due to Wernitz may be used. 

2. Traction improves with contacts whose major axis is in the rolling direction. 

3. Sideways forces due to the spin on the contact can be a substantial fraction of 
the transmitted traction. 

4. The effect of spin is always deleterious to traction efficiency. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 11, 1979, 

511-58. 
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APPENDIX A 


CORRECTIONS FOR DISK COMPLIANCE 

One of the implicit assumptions made in the present analysis is that the two con- 
tacting bodies are infinitely stiff in the plane of contact. This is not so in many practical 
situations. The effect of disk compliance may be studied by the following very simple 
model. Consider the layer of film between the two M contacting" disks. The combin- 
ation disk, film, disk form a sandwich of three deformable elements and relative ve- 
locities between the two disks is accommodated by shear in all three elements. When 
the film is thin and stiff, most of the shear may in fact take place in the disks. This is 
certainly the case when the disks roll together without a lubricant. For two bodies that 
roll together without a lubricant film it has been shown in references 12 and 13 that the 
following creep due to the transmitted traction force F occurs: 

(—) = q(-~ (ad 

VWdry V abG s/ 


where 

G„ shear modulus of disk material 

q constant depending on aspect ratio (given in appendix B). 

For an elastic lubricant film of thickness h under shear due to the transmitted traction 
force F the creep is given by 


AU \ = 32 h / F \ 
' U ' film 9rr a 


where it was assumed that 


G = 





(A2) 


( 8 ) 


The total creep for the sandwich under the imposed traction F 
of the disk creep and the film creep, namely, 


is given by the sum 


II 

I 

F 

(S L + 32hl\ 

V U /total 

7rab 

\ G g 97r a G/ 


(A3) 
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The parameter q covers the range of q = 0. 4 for k=0.2to q = 0. 93 for k = 5 
(longitudinal creep only). For commonly used traction drives operating at pressures 

_ O 

above 1. 2 gigapascals, h/a =* 10 and steel is generally employed as the disk material. 
When these values are used in equation (A3) together with a fluid shear moduli G of 
about 10 9 pascals, as found experimentally, then it is observed that most of the creep 
in the sandwich in the small slip region takes place in the disks, not in the fluid film. 
Reported shear moduli, as obtained from disk machine experiments have to be corrected 
for this effect (see, for example, refs. 5, 11, 14, and 15) for compliance corrected 
shear moduli. 

If, in the region of operation of traction drives, most of the elastic effects are 
caused by the disk material, then we are justified in using the shear modulus as a con- 
stant over the contact area since G a is not affected to any extent by the normal pres- 
sure. 
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APPENDIX B 


CORRECTIONS FOR CONTACT AREA, ASPECT RATIO, AND FILM THICKNESS 

Frequently traction data is available for a fluid at the same speed, temperature, 
and pressure but not at the same contact area because the roller geometry may have 
been different. From the analysis in the text it is clear that roller compliance and 
contact area play a significant role in the determination of the traction slope. An exact 
analysis for the correction of these terms is not yet possible; however, we can attempt 
the following simplified procedure. 

If we let m' be the slope of the traction curve for dry rolling bodies, then it fol- 
lows that the compliance corrected modulus for the fluid is approximately 

,bi> 


where 


G c compliance corrected modulus 
G apparent modulus 
m wet slope for contact 

m' dry slope for contact 

From the main text we have an expression for the apparent modulus equation 



Also from equation (Al), we have for the dry slope that 


F 

N 



7rabG 0 

o 

qN 


( 10 ) 


(B2) 


Now because N = (2/3) 7rabP Q we can simplify this equation a bit further to 


, 3 G s 

m’ = 


2 qP r 


(B3) 
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where q is aspect ratio dependent. Also equation (10) may be put in a slightly simpler 
form by this method 


G = 


P h 

Im-1 


4 a 


(B4) 


If we now assume that the compliance correction according to equation (Bl) is good to a 
fair degree, then for different measurements at the same pressure, speed, and temper- 
ature we might at least expect that the compliance corrected modulus would be the 
same; or if the * parameters refer to one set of experiments and the nonsuperscripted 
parameters to another set, then we may write under the aforementioned assumptions 
the following: 


G = g; 
c c 


By using equations (Bl) and (B4) this becomes 


Jl-SLY 1 

4 a \ m ' } 4 a \ m ' */ 


(B5) 


The comparison is made at the same pressure P Q so we can simplify the aforemen- 
tioned expression to 


*, * 

m h a 
mha* 

From equation (B3) for the same material and contact pressure we can write m'/m' = 
q* / q so that the final form of equation (B6) is 



m*h*a _ 
mha* 



(B7) 


Equation (B7) now relates the unknown slope m* to the other known parameters. 

We do need the parameters q (the Kalker coefficients) . These may be obtained 
from reference 13, and are plotted for both longitudinal traction measurements q x 
and side slip traction measurements q^. in figure 19. 

Equation (B7) is plotted in figure 20 for a common range of variables. 
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Figure 1. - Plan view of typical EHL contact indicating major 
variables; x is rolling direction. (The semimajor and minor 
axis are calculated from Hertzian theory.) 



Figure 2. - Sideview of typical EHL contact with constant separation 
between two contacting bodies. Pressure distribution in contact 
area is assumed to be Hertzian. 
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(a) Plan view of typical EHL contact showing inlet regions where 
shear stresses are assumed to be zero. 



distribution along y-axis in 
rolling direction due to spin 
and slip. Note location of 
speed pole e = AU la. 

Figure 3. - Plan view of typical EHL contact with its associated center line strain rate 
distribution due to slip, spin, and spin pole offset. 
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Traction stress, F v /7ra, MPa Dimensionless traction, F y /7TT 0 a' 



Slide slip/roll ratio, AV/U 


Figure 4 - Tests on hydratorque medium fluid in combined side slip and spin. Velocity in 
x-direction U = 1.9 meters per second; maximum Hertzian pressure in contact P 0 = 1. 3 
gigapascals; inlet temperature, 17° C. Elastic-plastic theory fitted to experimental results. 


Maximum 

Hertzian 


pressure, 



Figure 5. - Influence of average and pressure dependent shear properties G and t q on 
spin fraction prediction (data from ref. 5). Fluid, Turbo 33; velocity in x-direction U = 1. 2 
meters per second; inlet temperature, 3C P C. 
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Traction in x-direction, 



Figure 6. - Typical traction curve showing maximum traction 
coefficient \i and traction slope m. 



Aspect ratio of 
Hertzian contact. 
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Slip in x-direction, 

(b) Aspect ratio of Hertzian contact k - 4. 00. (c) Aspect ratio of Hertzian contact k = 8. 00. 

Figure 9. - Concluded. 






o 


e 

=3 

E 


.1 .2 .4 .6 .8 1 2 4 6 8 10 

Aspect ratio of Hertzian contact, k 

Figure 12. - Variation of maximum side force (J 5 ) with aspect ratio for 
contact under spin only. 



Speed pole, 


(a) Aspect ratio of Hertzian contact k = 1. 00. 

Figure 13. - Dimensionless spin torque J 6 versus the speed pole location J j/J 3 at various spin J 3 
val 
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(b) Aspect ratio of Hertzian contact k » 4. 00. 
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Slip and spin losses, (J 6 x Ji + J 6 x Jo} 


Region III - Rigid-plastic with spin 


Region II - Elastic-plastic with spin 


5 percent increase in at J 4 c a 75 compared to zero spin 

traction curves based on elastic-plastic analysis in fig. 7 

5 Percent increase in z at J 4 » 0.75 compared to rigid-plastic 

traction curves in fig. 8 


Region I - Elastic-plastic with no effect of spin 
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Figure 14. - Regions of influence of elastic and spin effects on slip Jj at constant traction value 
J 4 -0.75. 
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Slip in x-direction, J 1 

Figure 15. - Variation of J 7 = (J 4 x Ji + J 6 x J 3 ) with slip Ji for aspect ratio of Hertzian 
contact k « 1.00. 
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